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ABSTRACT 

In earlier work, we showed for the first time that the resolved stellar disk of NGC 300 is very extended with 
no evidence for truncation, a phenomenon that has since been observed in other disk galaxies. We revisit the 
outer disk of NGC 300 in order to determine the metallicity of the faint stellar population. With the GMOS 
camera at Gemini South, we reach 50% completeness at ig', i') = (26.8-27.4,26.1-27.0) in photometric 
conditions and 0.7" seeing. At these faint depths, careful consideration must be given to the background galaxy 
population. The mean colors of the outer disk stars fall within the spread of colors for the background galaxies, 
but the stellar density dominates the background galaxies by ^ 2 : 1 . The predominantly old stellar population 
in the outer disk exhibits a negative abundance gradient - as expected from models of galaxy evolution - out 
to about 10 kpc where the abundance trend changes sign. We present two scenarios to explain the flattening, or 
upturn, in the metallicity gradient of NGC 300 and discuss the implication this has for the broader picture of 
galaxy formation. 

Subject headings: galaxies: abundances — galaxies: individual (NGC 300) — galaxies: stellar content — 
galaxies: structure 



\ 1. INTRODUCTION 

Simulations of galaxy formation in cosmological context 
reveal that the process of galaxy assembly is expected to 
leave an imprint on the characteri stics of stars in the fain t 
outer regions of disk galaxies ( Bulloc k & Johns ton 20051). 
Due to their long dynamical timescales these regions retain 
the fossil record from the epoch of galaxy formation in the 
form of spatial and kinematic distributions, ages and chemical 
abundances of their stars (Freeman & Bland-Hawthorn 2002). 
These are also the regions where the effects of ste llar ra- 
dial mixing are expected to be mo st prominent (Roskar et al.l 
! I2008t ISchoenrich & Binnevll2008l) . The structure of the outer 
parts of galactic disks is therefore central to our understanding 
of the formation and evolution of disk galaxies. 

The studies of galactic disks have identified thre e basic 
classes of surface brightness profiles in galaxies (P ohlen et al.l 
l2007h . Truncations in surface brightness profiles, character- 
ized by the smooth break between the inner shallower and 
the outer steeper exponential, have been known for more 
than two de cades and detected in a large number of galaxies 
jvan der Kru it 1979; de Grijs et al. 2001; Pohlen et al. 2002). 
Despite this, the origin of the break remains unclear. The 
currently favored model for the explanation of these "sub- 
expon ential" truncations is the star formation threshold sce- 
nario (Kennicutt 1989'; 'Schavell2004l: iRoskar et al.l |2008', but 
see van der Kruit 2007). However, the threshold gas densities 
predi cted in these models are low (3- 10 M© pc~^; ISchavd 
120041) compared to the observed surface brightnesses at the 
radius of the break which co rrespond to stellar mass den- 
sities of 20-2000 Mq pc-2 dErwin et al.ll2007l) . The sec- 
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ond class of surface brightness profiles, so called "upbend- 
ing" or "super-exponential" breaks, have been explained by 
a scenari o involving mino r mergers, supported both observa- 
tionallv dlbata et al j[2005l) and by N-body/SPH simulations 
(lYoungeret al.ll2007l) . Until recently, it was believed that all 
galaxy disks undergo truncation near the Holmberg radius 
(26.5 mag arcsec"^) and that this was telling us something 
important abou t the c ollapsing protocloud of gas that formed 
the early disk (Ivan de r Kruit 1979) . This picture was ch al- 
lenged by our earlier finding (Bland -Hawthorn et al.ll2005h of 
a classic exponential disk (Freemanl ll970l) with no break in 
NGC 300, out to the distances corresponding to 10 disk scale 
lengths in this low-luminosity late type spiral. 

The existenc e of stars at such low stellar surface densities 
(0.01 Mq pc~^; lBland-Hawthorn et al.ll2005l) in the outer spi- 
ral disk has provided the motivation for this work. We seek to 
learn which physical process are responsible for the formation 
of the outer disks, and what the nature of the stellar popula- 
tions in this rarefied medium is. To this end, we obtained 
additional multi-band photometry (§',/') with Gemini/GMOS 
which allowed us to explore ages and metallicities of the stars 
in the outer regions of NGC 300. 

A number of studies have add ressed the question of s t ar for- 
mation in outer galactic disks. iDavidgd (I2003L l2006i l2007h 
find young and intermediate-age stars at large g alactocentric 
distances in NGC 2403, NGC 247 and M33. iBarker et al.1 
(12007) discovered young (less than 500 Myr old) main se- 
quences (MS) stars in the outer regions of M33, beyon d 
the break radius at 4.5 scale lengths dFerguson et al.ll20()7h . 
The existence of stars less than 0.5 Gyr old would sug- 
gest a very recent star formation episode in the outer disk 
of M33 as there probably has not been sufficient time for 
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these stars to be scattered from the inner disk and migrate 
outward. The same study also finds that the scale lengths 
of stellar populations in M33 increase with age. This is 
in contrast to the common scenario of inside-out galaxy 
formation in which the scale length of stars increase as 
the disks builds up to its present size, and it is expected 
that the mean age of the stellar disk decreases with radius 
jLarson' '1976"; 'Matteucci & Francois' '1989^, 'Chiappini et aJj 
[1997; Naab & Ostriker 2006; Muiioz-Mateos et al. 20oC 
These findings emphasize how little is understood about the 
formation of outer disks of spirals and call for deep photo- 
metric data on more pure disk galaxies. 

The Sculptor Group of galaxies, of which NGC 300 is a 
member, is one of the nearest and largest ensembles of disk 
galaxies outside the Local Group. Its galaxies are believed to 
be more or less isolated and as such are an ideal sample for 
our proposed study. NGC 300 is an almost pure disk galaxy 
(bulge light < 2%), with mild inclination (/ = 42°), which 
makes it well suited for stellar content studies. It lies at high 
galactic latitudes and therefore has low foreground reddening 
of E{B-V) = 0.01 1-0.014 mag (Schlegel et al. 1998), which 
translates into E(g'-i') = 0.018-0.024 mag. The high galac- 
tic latitude of Sculptor galaxies is also important because the 
contamination from the foreground Galactic stars is minimal 
at these high latitudes. In many respects NGC 300 is similar 
to the local spiral M33 (Blair & Long 1997). The o ptical disk 
scale length has been mea sured as 2 . 23 kp c i n B i (ICarignar 



f85l) and 1.47 kpc in I jKim et al.l l2004l). iTikhonov et al 
)05) and iButler et all (l2004 find that the stellar popula- 
tion is predominantly old and metal-poor, with no evidence 
for significant change in metallicity with time. The study of 
iPuche et all (Il990l) finds in NGC 300 an extended H I disk, 
which is severely warped in the outer parts (10' < r < 20'). 
The warp in H I d isk has often been l inked to the truncations 
in the stellar disk dvan der Kruitil2007h . 

We adopt a distance modulus of 26.5 1 for NGC 300. This 
is a weigh ted-mean of distances determined from the Cepheid 
Variables dFreedman et al.ll2001l; iGieren et al.ll2004l) and the 
tip of the red giant branch dButler et al.ll2004l; ITikhonov et al.l 
l2005h . 

The layout of the paper is as follows. We summarize the 
details of the observations and the data reduction in ^ and 
describe the photometric measurements and artificial stars test 
in ^ The color-magnitude diagrams, star counts and surface 
brightness profiles, and the metallicity distribution function of 
NGC 300 are presented in ^ Finally, we discuss our results 
inS 

2. OBSERVATIONS AND DATA REDUCTION 

The deep g' and /' images of three fields in the outskirts 
of NGC 300 were obtained using the Gemini Multi Object 
Spectrograph (GMOS) on Gemini South telescope as a part 
of the program GS-2005B-Q-4; in total 34 hours of data have 
been taken. The observations are summarized in Table [T] The 
NGC 300 fields span a major axis distances of 7- 16 kpc 
and were made to overlap slightly to allow for the check of 
the relative photometric calibration between the fields. The 
inner-most field was placed at the optical edge identified in 
the Digitized Sky Survey (DSS) to provide continuity with 
earlier work. The locations of the fields are marked in the 
Figure[T] The field of view in a single GMOS image is 5'.5 on 
a side; at the distance of NGC 300 1' corresponds to 0.6 kpc. 

The GMOS detector is a mosaic of three 2048 x 4608 LEV 
CCDs dHook et alj|2004l) . The raw pixels were binned 2x2 
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Fig. 1.— The Digitized Sky Survey (DSS) wide-field image of NGC 300. 
(Nortli is to the top, east is to the left.) Squares mark the positions of GMOS 
5'.5x5'.5 field of view. 

during readout to provide a better match of the image sam- 
pling to the seeing, producing the pixels of 0". 146 in size. 
The average image quality of the data is 0".7 FWHM in g' 
and 0".6 FWHM in /'. 

For each field, 13 x 600 s exposures were taken in g' and 
22 X 600 s exposures in /' band. Our goal was to achieve 
3(7 photometry in g' and /' for a G giant star with (V= 27.5, 
B= 28.5, 1= 26.5). The data were obtained in a dithered 
pattern, which allowed for the gaps between the individual 
GMOS CCDs to be filled when constructing the final image. 

The data were reduced using the standard IRAF/Gemini 
routines. Due to the known problems with GMOS bias 
frames, we inspected all bias frames and discarded those in 
which the counts varied significantly between the subsequent 
frames. For the same reason, we decided to have the over- 
scan regions subtracted in the reduction process. The process- 
ing pipeline included creating the master bias and flat field 
frames (gbias, gif lat), bias subtraction and flat fielding 
(gi reduce), mosaicking of individual GMOS CCDs into a 
single reference frame (gmosaic), and finally combining the 
dithered exposures into a final image (imcoadd). In addi- 
tion, /' data were affected by fringing and the fringes needed 
to be removed. The master fringe frame was created for each 
field from the individual reduced image frames (gif ringe) 
and subtracted from the individual images before coadding 
(girmf ringe). 

After the reduction procedure, the artifacts of mosaicking 
the GMOS CCDs into a single image were still present in g' 
data, visible as a background offset between the chips in the 
final images. The level of the offset was of the order of ~ 
1% of the background, but was clearly evident upon visual 
inspection, as shown in the left panel of Figure |2l In order to 
attempt to remove these effects, we applied the same fringe 
removal procedure on our g' data as described above for /' 
images. The procedure improved the background consistency 
between the chips; the example of the result for one of the 
fields is shown in the right panel of Figure |2l 

3. PHOTOMETRY AND COMPLETENESS ANALYSIS 

3.1. Photometry 

The final g' and /' band GMOS images we re analyzed us - 
ing die DAOPHOT/ALLSTAR software suite dStetsonlll987l) . 
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TABLE 1 
GMOS Observing Log. 



Target 


RA (J2000) 


DEC (J2000) 


Date 


Filter 


Exposure (s) 


FWHM (") 


NGC300 1 


55 50 


-37 47 18 


2005 Sep 7-8 


g' 


13x600 


0.71 




55 50 


-37 47 24 


2005 Aug 4 


i' 


22x600 


0.66 


NGC3()0 2 


56 15 


-37 46 45 


2005 Aug 8-10 


g' 


13x600 


0.72 




56 15 


-37 46 45 


2005 Aug 8-10 


i' 


22x600 


0.67 


NGC300 3 


55 59 


-37 50 59 


2005 Oct 10-11, 31 


g' 


13x600 


0.82 




55 59 


-37 50 53 


2005 Sep 8, Oct 8-10 


i' 


22x600 


0.72 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and 
arcseconds. 




Fig. 2. — The reduced g' images of the NGC 300 F2 field. Left panel: the image after the standard reduction procedure has been applied (see text for details). 
Right panel: the image after applying the additional fringe removal procedure. This additional reduction step significantly improved the consistency of the 
background between the three GMOS CCD chips. 



For each field DAOPHOT produced object catalogues and 
performed aperture photometry in a series of circular aper- 
tures. About 100-200 relatively bright and isolated stars 
in each field were chosen as PSF stars and used to itera- 
tively compute the point spread function (PSF) for the field. 
All PSF stars candidates were examined visually within the 
IRAF/daophot package and their radial, contour and mesh 
profiles were inspected for any indications of their non-stellar 
nature. After the PSF stars were subtracted within the stan- 
dalone DAOPHOT suite, the positions of subtracted PSF stars 
were inspected again using the above indicators. All PSF star 
which were not subtracted cleanly were excluded from the 
subsequent PSF calculations. In addition, stars which had PSF 
subtraction errors more than 3a away from the mean value, 
as reported by standalone DAOPHOT package, were also ex- 
cluded. The next iteration PSF was calculated from the image 
in which, within the fitting radius of each PSF star, all but 
PSF stars were subtracted. The procedure described above 
was repeated once more to derive the PSF; finally, ALLSTAR 
was used to fit this PSF to the objects in the catalogues and 
determine their PSF photometry. 

In addition to the information on the object magnitude and 
the related error, ALLSTAR provides information on the qual- 
ity of PSF fitting for each star through the parameters CHI 
and SHARP. CHI represents the ratio of the observed pixel- 
to-pixel scatter in the fitting residuals to the expected scatter. 
It is expected for the values of CHI to scatter around unity 



when plotted against t he derived magnitude (cf. Figure 28 of 
IStetson & Harts|[T988 ). In our analysis however, the values 
of CHI are scattered around 0.4. This may be indicative 
of the incorrect values for the gain and readout noise used in 
the analysis (P. Stetson, private communication). The corre- 
spondence with Gemini staff confirmed that the gain and read- 
out noise values used are indeed correct and the reason for 
the invalid CHI values in ALLSTAR remained unclear. How- 
ever, any photometric errors introduced by potentially incor- 
rect gain and readout noise values will be taken into account 
through our artificial star tests which we use to estimate the 
observational uncertainties in the data. 

The aperture corr ections were calculated using the 
DAOGROW software (IStetsonlll99^ . which performs a so- 
phisticated growth-curve analysis to derive a "total" magni- 
tude for each star. For each field, multiple-aperture photome- 
try was done on the image frames from which all but PSF stars 
have been subtracted. The weighted mean difference between 
the PSF-based magnitude (from ALLSTAR) and the "total" 
magnitude (from DAOGROW) of the PSF stars was adopted 
as "aperture correction" and this correction was applied to 
PSF-based magnitudes of all stars in the frame. The uncer- 
tainties in the aperture correction over all bands and fields 
were 0.006-0.008 (the quoted uncertainty is the standard er- 
ror of the mean). 

The photometric calibration for the GMOS data was es- 
tablished using the iLandolj ( Il992h standard star fields. The 
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Fig. 3. — Completeness curves for the observed NGC 300 fields, in g' 
(top panel) and i' (bottom panel). The horizontal dashed line marks the 50% 
completeness hmit. 

g' and /' magnitudes of standard stars were calculated using 
the relations from Fukugita et al. ( 1996). In total, 6 Landolt 
fields were observed multiple times over the two months run, 
with the total of ^ 140 standard stars observations. The total 
observed magnitudes of standard stars were measured using 
DAOGROW and the zero point was calculated as the differ- 
ence between the observed magnitude corrected for the expo- 
sure time and the standard magnitude. Zero points varied by 
0.02 between the standards stars observations on the different 
nights of the same month, or by 0.12 and 0.07 in g' and /', re- 
spectively, over the two month period. Photometry was trans- 
formed onto a single calibration by comparing the magnitudes 
of stars in the overlapping regions, taking as a reference the 
field with the best calibration of the three (Fl in i' and F2 in 
g', both observed during the August 2005 run). Results of the 
completeness analysis described below were used to confirm 
that the scatter of the photometry in the overlapping regions 
agrees with the scatter predicted from the artificial star tests. 

3.2. Completeness Analysis 

When interpreting stellar photometry of faint stars in 
crowded fields, it is important to understand the observational 
uncertainties inherent in the data. It is a standard procedure 
to use artificial star tests to estimate the completeness of the 
data and assess the effects the crowding has on the calculated 
magnitudes and their errors. In the artificial star test a set of 
stars with known magnitudes is added to the frame, and the 
frame is analyzed using the standard data reduction pipeline 
to recover the photometry of all stars, real and artificial. 

To be able to estimate the data completeness robustly, a 
large number of stars has to be added to the images and an- 
alyzed. This number should be significantly larger than the 
number of real stars detected in the image. At the same time, 
it is important not to affect the crowding of the images signif- 
icantly, as that will result in the unrealistic estimates of com- 
pleteness. As a compromise between the crowding in the field 
(where a small number of artificial stars is optimal) and the 
computing efficiency (where a large number of artificial stars 



is favorable as it decreases the total number of images that 
need to be analyzed), the number of stars that we add to each 
frame is 5 - 10% of the number of detected stars in the frame. 
We then run a hundred tests on each image to obtain statis- 
tically robust estimate of the completeness and observational 
uncertainties. 

In order to sample the crowding conditions in the frame cor- 
rectly, artificial stars injected into the image must not overlap. 
To make certain that there is no overlap between the artificial 
stars, we add them to the frames in a regular grid. The grid 
has 30 stars on a side (i.e. 900 in total), with a separation of 
69 pixels (10") between neighboring stars. Our choice of the 
grid and cell size is guided by the considerations explained 
above and the need for the whole extent of the image to be 
sampled in each run. A total number of stars added to each 
frame over 100 artifical star test runs is ^ 10000 per frame. 
The grid origin has been randomly offset between the runs in 
order to enable sampling of the whole extent of the image. 
The same grid setup is used to add artificial stars to one (§',/') 
pair of images. 

To assess the color effects it is important that the colors 
assigned to artificial stars correspond to the realistic stellar 
populations. To ensure that the magnitude and color distri- 
bution of our artificial stars matches magnitudes and colors 
in the original data frames, we construct artificial photome- 
try by randomly selecting stars from the original CMDs and 
assigning their g' and /' magnitudes to artificial stars. 

Artificial stars with the positions and magnitudes deter- 
mined as described above are added to the frames using the 
DAOPHOT/ADDSTAR routine. The resulting images are an- 
alyzed with the same data reduction pipeline as the origi- 
nal science frames, using the PSF calculated for each field 
as described above. A star is considered recovered if it is 
detected in both g' and frame and if the difference be- 
tween the injected and recovered magnitude in both bands is 
< 0.5 mag. Completeness as a function of magnitude is cal- 
culated as the ratio of the number of recovered and injected 
stars in a given 0.5 magnitude bin. Completeness curves for 
three NGC 300 fields are shown in Figure |3] We compute 
the completeness in terms of the instrumental magnitude and 
then transform them using the appropriate photometric cali- 
bration equations into the completeness in terms of the stan- 
dard magnitude. Our photometry is 50% complete down to 
(§ ',/')=(26 .8 - 27 .4,26 . 1 - 27 .0). 

To check the accuracy of the photometry we calculate the 
difference between magnitudes of injected and recovered ar- 
tificial stars. Figure in shows the difference between the input 
and output magnitudes as a function of the input magnitude 
for one of the observed fields. Over the most of the magni- 
tude range in g' and /' the photometric error with which the 
stars are recovered is < 0.3 mag. There is a systematic trend 
over the whole magnitude range, with stars being recovered 
brighter than their true magnitude. Assuming that this effect 
is purely due to blending, we estimate the blending fractions 
to be 5- 10% in F2 and F3 and ~ 20% in Fl . However, the cal- 
culated blending fractions are an overestimate because at the 
faint end random fluctuations in the unresolved background 
may scatter faint stars below the detection limit, with only the 
stars lying on the positive background fluctuations being re- 
covered. 

We use artificial star tests to derive uncertainties of our pho- 
tometry. We adopt as photometric errors the mean of the abso- 
lute value of the difference between the input and output mag- 
nitudes of artificial stars. Figure |5] shows the errors reported 
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Fig. 4. — The difference between the injected and recovered magnitudes of 
all recovered artificial stars, in g' (top panel) and i' (bottom panel), for one 
of the observed NGC 300 fields. Only the stars with the difference between 
the input and output magnitude smaller than 0.5 were considered recovered 
when calculating completeness curves shown in Figure |3] Other NGC 300 
fields show similar behavior 
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Fig. 5. — Photometric errors as given by ALLSTAR (grey points) and 
calculated from the artificial star tests (black squares); top panel shows the 
results for g' band, bottom panel for band data. Photometric uncertainties 
are derived from the artificial star tests as the mean difference between the 
input and output magnitudes of artificial stars in 0.5 magnitude bins. Error 
bars give the scatter in each bin. 

by ALLSTAR (points) and photometric uncertainties calcu- 
lated from the artificial star tests as described above (squares), 
as a function of input magnitude. The error bars are the scatter 
in each 0.5 mag bin. At the faint end the errors derived from 
the artificial star tests are ^ 0.04 mag larger than the errors 
given by ALLSTAR. 

In addition to estimating completeness and photometric er- 
rors, we use the artificial stars tests to discard the non-stellar 
and spurious detections. Due t o the non-standard CHI output 
of the ALLSTAR routine (see ^3.1b . we are not able to apply 
commonly used cuts in CHI parameter to reject false detec- 
tions. With the artificial star tests however, we are able to use 
the information on CHI and SHARP of the artificial stars to set 
the criteria for removing the false detections from the data. In 
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Fig. 6. — Values of the CHI (bottom panel) and SHARP (top panel) param- 
eters as reported by ALLSTAR for all recovered artificial stars (grey) and the 
real stars (black). Dashed lines are the cuts in CHI and SHARP used to discard 
spurious detections. 

Figure|6]we show the distribution of CHI and SHARP parame- 
ters for the real stars (black points) and all recovered artificial 
stars (grey points) as a function of magnitude, for one of the 
observed NGC 300 fields. In order to remove non-stellar ob- 
jects from the sample, we employ cutoffs based on the values 
of SHARP parameter of the recovered artificial stars. We di- 
vide the data into 0.5 magnitude wide bins and discard objects 
that are more than 3a away from the mean value of SHARP in 
a given bin. The dashed curve in the top panel of the Figure|6] 
shows the extent of the allowed SHARP values. In addition, 
we use a constant CHI limit and remove from the stellar cata- 
logues all stars with CHI > 0.5. 

4. RESULTS 
4. 1 . Color-Magnitude Diagram 

We show the {g' - i',i') color-magnitude diagrams (CMDs) 
of the three observed NGC 300 fields in Figure |7] The most 
prominent component of these CMDs is the well-populated 
red giant branch (RGB). This feature is typical of a system 
composed of a stellar population older than 1 Gyr. In all ob- 
served fields the red giants branch stars make up ~ 50% of 
all stars detected in the field. A small number of red stars is 
found above the tips of the theoretical red giant branch tracks 
(Figure HJ; these stars could be intermediate-age asymptotic 
giant branch stars and young red supergiants. Also marked on 
the CMDs in Figure|2]are the 50% completeness limit and the 
photometric errors at a given /' magnitude and g' - i' color of 
1, as determined from the artificial stars tests 03.2b . 

In addition to the dominant old and red stellar population, 
we detect a sparsely populated main sequence with g' - i' col- 
ors from -1 to ~ 0. This feature contains 10% of the stars 
detected in the outer disk of NGC 300 and most likely rep- 
resents main-sequence stars with the ages of a few hundred 
Myr 

4.2. Star Counts Profile 

In ou r previo us study of the outer regions of NGC 300 
(iBland -Hawthor net al.ll2005h we used r' star counts to trace 
the faint stellar disk out to 15 kpc (25')- Now with the 
multi-band photometry in hand we can refine our radial pro- 
file to target distinct stellar populations in the outer disk of 
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Fig. 7. — Color-magnitude diagrams (CMD) of the three NGC 300 fields. EiTor bars show the photometric uncertainties from the artificial star tests at a given 
i' magnitude and the color of g' — i' = I. The dashed lines indicate the 50% completeness limit. 
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Fig. 8. — Colo r magnitude diagram of the three NGC 300 fields, with the 
Victoria-Regina I'Vanden Berg et alj |200fll) set of [o/Fe]= 0.0, 8 Gyr stellar 
evolutionary tracks overplotted. The model grid has the metalHcities in the 
range [Fe/H]=[-2.31,0.00]. Two regions marked in the figure were used to 
select stars belonging to red giant branch (RGB) and main sequence (MS) 
population. 

NGC 300. 

We derive the radial profile of the red giant branch and main 
sequence stars by counting the stars in elliptic annuli with 
the ellipticity (inclination) and position angle corresponding 
to that of NGC 300 (PA= 1 10°, / = 42°). In Figure[8]we repro- 
duce the color-magnitude diagram of the NGC 300 Fl field 
and define the selection boxes we use to select RGB and MS 



stars. To account for the crowding effects we compute com- 
pleteness as a function of galactocentric distance using the 
output of the artificial stars tests described in 33.21 This is cal- 
culated as a ratio of the number of recovered to injected stars 
in each 0.5' radial bin and used to correct star counts pro- 
files. The resulting completeness-corrected profiles of RGB 
and MS stars are shown in the top panels of the Figure|9] 

In order to be able to derive the intrinsic star counts profile 
it is crucial to as reliably as possible estimate the contribution 
of the contaminant sources - dwarf stars in the Galaxy and 
faint unresolved background galaxies that get mistaken for 
NGC 300 stars. Given the high galactic latitude of NGC 300 
(b = -79.4°) we do not expect the contribution from galactic 
dwarfs to be significant. Besancon models (Robin et al. 20031) 
predict 26 stars per GMOS field of view in the magnitude and 
color range corresponding to our RGB selection box; this ac- 
counts for 0.5 - 1 % of all detected objects in a single field with 
the same color-magnitude criteria. 

The online Gal axyCount calculator 

(Ellis & Bland-HawthornI '2007) is a useful tool for esti- 
mating galaxy counts in deep wide-field images. It had 
confirmed our earlier finding that even when the contribution 
from the faint unresolved galaxies is accounted for, the 
exponential disk of NGC 300 as traced by r'-band star counts 
extends out to 1 disk scale lengths dBland-Hawthorn et al J 
,2005; Ellis & Bland-Hawthoriil l2007h . However, we now 
derive star counts profiles for different stellar populations 
which populate distinctive regions in the color-magnitude 
space (Figure [8]) and require information on the color, as well 
as magnitude, of the background galaxy population. Since 
GalaxyCount at present does not include the color infor- 
mation, we can only estimate the number counts predicted 
in g' and /' bands separately. For example, in the case of 
our RGB selection box (Figure [U GalaxyCount predicts 
82- 114 galaxies in the magnitude range g' = 23.9-28.0 and 
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Fig. 9. — Star counts profile of the red giant branch (left panels) and main sequence stars (right panels) in the outer disk of NGC 300. Stars inside the regions 
marked in Figure[8]were used to calculate the profile. The top panels show the profiles before the contribution of the faint unresolved galaxies has been subtracted. 
The dashed line corresponds to the predicted number of background galaxies (per arcmin^) with the colors and magnitudes within our RGB and MS selection 
boxes. The bottom panels show the star counts profiles corrected for the background galaxies contribution. 

114- 140 galaxies in the range i' = 22.7-26.4. To estimate 
the galaxy number counts we use the completeness curves 
derived separately for each field; the range in the number 
counts estimates reflects the differences in the parameters 
describing the completeness curves for different fields. 
However, the number counts predicted by GalaxyCount 
are an overestimate of the true contamination within our RGB 
selection region because (1) they are derived without any 
restrictions on the color of the background galaxy population 
and (2) a number of contaminating galaxies has already been 
somewhat reduced by the use of the SHARP statistic. 

For the reasons described above we chose to use the 
galaxy counts from the William Herschel Deep Field (WHDF, 
iMetcalfe et al. 2001) to estimate the contamination from the 
faint background galaxy population. The main caveat with 
this kind of approach is that the WHDF survey is not as deep 
as our data and some assumptions need to be made in order 
to estimate the number counts of galaxies in our data using 
a shallower survey. We calculate the /'-band number counts 
of the galaxies falling within our color selection criteria and 
fit linearly the (log of) differential number counts in 0.5 mag 
bins in the range covering 3 magnitudes above the WHDF 
magnitude limit. To determine the galaxy counts below the 
limit of the WHDF survey we assume that the counts in the 
bins 2-3 magnitudes below the survey limit follow the same 
linear trend (in the log space) as the counts in the brighter bins 
used in the fit. We finally take into account the incomplete- 
ness of our own data and correct the derived galaxy number 
counts using the completeness curves presented in Figure [3] 
The estimated mean number of unresolved background galax- 
ies over all fields is 36 ± 7 per arcmin"^ for our RGB selection 
box and 15 ± 8 per arcmin"^ for the region used to select the 
MS stars. The errors reflect both the variation of the estimated 



galaxy counts between the three fields and the cosmic vari- 
ance. The latter was estimated using GalaxyCount, which 
for given galaxy number counts in a given magnitude range 
calculates the expected variance. The star counts profiles be- 
fore and after subtraction of the background galaxy popula- 
tion are shown in Figure |9] in top and bottom panels, respec- 
tively. Error bars in the bottom panel are a combination of 
errors quoted above and Poisson uncertainties. Under the as- 
sumption that the slope of the WHDF number counts changes 
at the survey magnitude limit to a value twice that derived 
from the bins 3 magnitudes below the survey limit, we derive 
the galaxy counts as 80 ± 15 and 29 ± 14 per arcmin"^ for our 
RGB and MS selection, respectively. After the subtraction of 
the newly determined galaxy counts, the old stellar disk re- 
mains purely exponential but falls below the detection limit at 
^ 20 arcmin. However, these higher counts are very likely to 
be an unrealistically high estimate for the contribution from 
the background galaxy population, given that they are signif- 
icantly higher than our MS star counts in the outermost two 
thirds of our MS radial profile. 

As the Figure|9]clearly shows, while the old stellar disk ex- 
tends out to 25', corresponding to 15 kpc, the young main se- 
quence stellar population is more centrally concentrated and is 
only detected out to 17' (10 kpc) along the southern end of the 
major axis. The disk, as traced by the old stellar population, 
does not exhibit a break in surface brightness profile once the 
star counts have been corrected for the contamination by unre- 
solved background galaxies. This complements our finding of 
an unbroken exponential disk in late type galaxy NGC 300 out 
to the distances correspondin g to about 10 disk scale lengths 
jBland-Hawdiorn et alj|2005l) . 

4.3. Surface Brightness Profile 
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Fig. 10. — Surface brightness profile of NGC 300. Open circles are I- 
band surface photometry data from Kim et'al] 120041) : the data points have 
been shifted downward by 1.2 mag. Full squares are Bj surface photometry 
measurements taken from Carignan 1 1985) and have been shifted upward by 
0.4 mag. The star symbols are our surface brightness measurements derived 
from star counts. 

In Figure [To] we derive the surface brightness profile of the 
outer disk of NGC 300 from the /'-band starcounts of the 
old RGB stars. The inner data points in this figure are taken 
from the s urface photometry studies by Kim et alj (12004) and 
ICarignanI {1985). The outer disk surface brightness is cal- 
culated in 0.5' annuli by summing the flux of all stars in 
a given annulus and correcting for the radial completeness, 
number of pixels in the annulus and the pixel scale to ob- 
tain surface brightness in the units of mag arcsec"^. We cor- 
rect for the background galaxy contamination by randomly 
removing from each radial bin a number of objects corre- 
sponding to our estimate of the galaxy number counts (36 
per arcmin"^) corrected for the area of the annuli. We repeat 
this process by drawing 100 random "background galaxy- 
corrected" samples from our full stellar catalogue and use this 
to calculate the final value of surface brightness in each ra- 
dial bin, and estimate its error bars (as standard deviations of 
the full set of 100 samples). Finally, we correct for the light 
missed due to the fact that we detect only the brightest stars 
in the outskirts of NGC 300 and the inclination of the galaxy. 
In iBland-Hawthorn et al.l (l2005h we estimate the fraction of 
missing fight to be ~ 45% for a study of a similar photomet- 
ric depth. This translates to intrinsic surface brightness 0.65 
mag brighter than what we derive here. The galaxy inclination 
(/ = 42°) somewhat counteracts this effect; if the galaxy was 
viewed face-on the resulting surface brightness would be 0.32 
mag fainter. Therefore, to account for both effects the cal- 
culated surface brightness has to be shifted upward by 0.33 
mag. 

The surface brightness profile in Figure [TO] emphasizes the 
power of resolved star counts over the traditional diffuse light 
imaging for reaching down to very low surface brightness 
levels necessary for studying the extremely faint outskirts of 
galaxy disks. 

4.4. Metallicity Distribution Function 

Given the color-magnitude diagram and the distance to 
NGC 300, we are able to compute the metallicity distribu- 
tion function of the stellar populations in its outskirts. The- 
oretical isochrones show that metallicity (rather than age) 
is the primary factor affecting the color of the red giant 
branch (VandenBerg et al. 2006; see also Figure 5 and 6 of 
[Harris et al. 1999); for example, at absolute magnitude M/ = 
-3 the 0.2 mag blueward shift in V - / can be achieved by 
a decrease in age from 1 3 to 6 Gyr, or by a relatively small 
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Fig. 1 1. — Metallicity distribution function of NGC 300 stars calculated 
using the Victoria-Regina stellar evolutionary tracks with the age of 8 Gyr 
and no a-enhancement. Also shown are the metallicity distribution functions 
calculated assuming the stellar ages of 12 Gyr and [a/Fe\= 0.0 (dashed line) 
and an 8 Gyr old stellar population with [a/Fe]= 0.3 (dotted line). If the 
stellar population in the outskirts of NGC 300 is a-enhanced or older than 
the assumed age of 8 Gyr, our procedure will have sUghtly overestimated the 
number of stars at the high-metallicity end, but without a significant effect on 
the mean metallicity of the population. 

decrease in [Fe/H] from -1 .0 to -1 .2 (iMouhcine et alJl2005h . 
This high sensitivity of photometric properties on metallicity 
allows us to derive the metallicity distribution function from 
the photometry of the red giant branch stars. 

To derive the metallicity distribution function we convert 
the observed distance- and reddening-corrected stellar pho- 
tometry to metallicity on a star-by-star basis. We superim- 
pose on observed color-magnitude d i agram the stellar evo- 
lutionary tracks of I VandenBerg et akl (|2006|) and interpolate 
between them to derive an estimate of star's metallicity (Fig- 
ure [8]|. The model grid consist of 16 finely spaced red giant 
tracks for a 8 Gyr stellar population without a-enhancement, 
covering the range of metallicities from [Fe/H]= -2.3 1 to 0.00 
in the steps of approximately 0.1 dex. We define a region in 
ig' -i',i') which we expect to be populated primarily by red 
giant branch stars and only consider the stars within this se- 
lection box when calculating metallicities (Figure[8]). 

In computing the metallicity distribution function we 
choose to use the theoretical stellar evolutionary tracks rather 
than empirical globul ar cluster fiducials. The g rid of red gi- 
ant branch tracks from lVandenBerg et al.l (l2006l) on which we 
base our metallicity estimates is much more finely spaced 
than often used Milky Way globular cluster fiducials by 
Saviane et all (120 00); as shown by iHarris et^ fl999 ) and 
Harris & Harrisl " i200Q) the coarseness of the model grid may 
introduce artifacts in the calculated metallicity distribution 
function. By inspection of the Figure[8]it is clear that the stars 
in the outer regions of NGC 300 are metal poor and there is 
no need f or adding stellar tracks m ore metal rich than those 
present in lVandenBerg et aTl (l2006h grid. 

In order to take into account the incompleteness effects we 
compute the completeness-corrected metallicity distribution 
function by counting each star (Cj'C,') ' times, where Cg' and 
C,' are completeness in g' and /' at a given point in the CMD. 
We only consider stars which fall above 50% completeness 
limit in both g' and i' . For fainter stars the completeness cor- 
rection factor becomes too large and not necessarily a realistic 
representation of the incompleteness effects which are poorly 
understood at these faint magnitudes. 

The presence of stars of different ages will affect our metal- 
licity estimates which utilize the 8 Gyr stellar tracks; the ef- 
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Fig. 12. — Metallicity as a function of radius for all stars in the outer disk 
of NGC 300. Full circles are the median metallicities in 1' bins. Solid Hues 
represent the metallicity gradient and are a linear fit to the median metallic- 
ities in the range R < 15.5' and R > 15.5'. Error bars represent the mean 
uncertainty in [Fe/H] in each bin. 

feet is illustrated in Figure [TT] The solid line is the metallic- 
ity distribution function calculated using the stellar tracks for 
the 8 Gyr stars with no a-enhancement. Also shown are the 
metallicity distributions constructed assuming an older age 
(12 Gyr, dashed line) and different value for a-enhancement 
([a/Fe]= 0.3, dotted line). If the true stellar population is 
older or more a-enhanced than our adopted stellar evolution- 
ary tracks, this procedure will have overestimated the num- 
ber of high-metallicity stars, and similarly, underestimated the 
contribution of low-metallicity stars, but the effect on the peak 
metallicity of the stellar population will not be significant. 

4.4.1. Metallicity Gradient 

In Figure [12] we plot metallicity as a function of depro- 
jected galactocentric radius for all stars in the outer disk of 
NGC 300. Preliminary analysis of stellar metallicities has 
shown that metallicity decreases with radius out to about 15.5' 
after which the gradient becomes mildly positive. We there- 
fore divide the data into R < 15.5' and R > 15.5' (9.3 kpc) 
bins and independently calculate the metallicity gradient (in 
dex kpc"') for each region, by fitting linear functions to me- 
dian metallicities in 1' wide radial bins: 

/;< 15.5' : [Fe/H] =-0.11(±0.02)/;-0.04( ±0.15) 

R > 15.5' : [Fe/H] =0.034(±0.004)/;-1.37(±0.06) 

Despite the large scatter in calculated metallicities, we 
clearly detect both the negative abundance gradient out to 
^ 9.3 kpc and a change in the trend outwards of 9.3 kpc. 

Figure [13] shows the pure exponential surface brightness 
profile of NGC 300 and its metallicity gradient on the same 
radial scale. In the bottom panel of this figure we compile 
the available data on the inner disk metallicities from the lit- 
erature. All metallicities except the recent Kudritz ki et al.l 
(i2008|) results refer t o oxy gen metalhcities on l y, and all ex- 
ceptl Kudritzkietal.1 (l2008h and lUrbaneia et al.1 ( 120051) metal- 
licities are derived using the H II -region emission lines. 

The reasons for the metallicities we derive being seemingly 
too low compared to the spectroscopic metallicities in the cen- 
tral regions of NGC 300 are two-fold. Firstly, the inner disk 
metallicities in the Figure [13] either represent gas, rather than 
stellar, metallicities, or were measured from young blue su- 
pergiants, both of which are expected to have higher metal- 
licities than the older red giant branch stars at the same radii. 
Secondly, as can be seen from Figure[T2] the bulk of the stars 
populate the high-metallicity end of the plot and it is the low- 
metallicity tail of the abundance distribution function (see 
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Fig. 13.— Top: Surface brightness profile of NGC 300. Bottom: Metal- 
licity gradient in the disk of NGC 300. Inner disk data points are from the 
spectroscopic abundance studies listed in the figure. Filled circles are outer 
disk binned mean meta llici ties and soHd lines are the linear fits to these values 
as calculated in Figure [12] 

Figure [TTT i that is responsible for lower mean metallicities. 
Since we base our discussion primarily on the overall shape 
of the metallicity gradient, rather than the absolute value of 
metallicities this effect will not influence our conclusions. 

In order to test the robustness of the derived abundance gra- 
dient, we performed a series of tests in which we randomly 
removed from the sample (i) a constant number of stars, cor- 
rected for the different areas of elliptic annuli, or (ii) a number 
of stars proportional to the number density of stars in a given 
radial annuli. These procedures mimic the effect of correcting 
for the (i) constant contamination (e.g. background galaxies, 
stellar halo) and (ii) radially dependent contamination (e.g. 
thick disk). Only after artificially removing 50% of stars in 
each bin are we able to derive a flat abundance gradient over 
the full radial extent sampled by our data. However, there is 
no realistic mechanism that could account for a 50% contam- 
ination level within our RGB selection region and hence we 
consider our derived abundance gradient to be robust. 

5. DISCUSSION 

In this paper we revisit NGC 300 to investigate the faint 
stellar population in its outer d isk. We confirm the earlier 
finding dBland-Hawthorn et al] r2005h of the extended stel- 
lar disk in NGC 300 with the surface brightness declining 
exponentially out to ^ 25', or about 10 disk scale lengths 
04.2I ). Only three other disk so far have been traced down 
in resolved stars to the same effective surface brightness lev- 
els- that of the Galaxy (ll bata et al. 2003), M31 ( Irwin et aP 
l2005h and M33 dFergusonet al.. .2007) . Except for the case 
of a "sub-exponential" disk of M33 the other three disks are 
either pure exponentials or "super-exponential" and are de- 
tected down to the effective surface brightness levels below 
r^ 30 mag arcsec"^. The existence of stars at these extremely 
low surface densities in the outer galactic disks is puzzling. 
We want to learn what the nature of the faint stellar population 
in NGC 300 is and whether its outer disk is young or old. The 
answers to these questions will have important implications 
for the scenarios for the formation of spiral galaxies. With 
the multi-band imaging in hand, we are now able to assess the 
ages and metallicities of stars in the outskirts of this late type 
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spiral. 

From the (g ' color magnitude diagrams of NGC 300 

presented in ^4.11 and the star counts profiles of RGB and 
MS stars (Figure |9]l we find that its outer disk is predom- 
inantly old. Old outer disks are in disagreement with the 
inside-out picture of spiral galaxy formation in which the 
inner disk is built up on shorter timescales than the outer 
disk and yo ung stars are expe cted to inhabi t the larger radii 
jLarsonl [l97 6; Matteucci & F rancois! 119891; IChiappini et alj 
119971: iNaab&O striker 2006). This is however not the first 
time the old stellar popula tion has been discove r ed at large 
galactocentric distances. iFerguson & JohnsonI (1200 Ih use 
WFPC2 on board HST to observe a field on the major axis 
of M31 at the radius corresponding to 3 disk scale lengths. 
Using the mean 1-band magnitude of the red clump as an 
age indicator they find that the mean age of the stellar pop- 
ulation in the outskirts of M31 is > 8 Gyr. In their stud- 
ies of the ou t ermos t regions of M33. lGalleti etall jlool and 
[Barker et al.l (l2007h both find that the young main sequence 
stars appear to be more centrally concentrated and hence have 
smaller scale lengths than the older red giant branch popu- 
lation. The disks that have formed inside-out however have 
stellar populations which get progressively younger with ra- 
dius and the oldest disk stars have the s mallest scale lengths, 
in contrast to what is observed in M33 (Matteucci & Francois 
1989; Chiappini et al. 1997; Galleti et al. 2004; Barker et al. 
20071). 

It may be poss ible to reconcile the two results in the con- 
text of models bv lRoskar et al.l (120081) . In their simulations of 
spiral galaxy formation the old outer disks arise as a conse- 
quence of stellar migrations from the inner disk. Star forma- 
tion rate density falls off sharply at the radius of a break in 
the radial profile, but the stars are found in the outer disk at 
least ^ 5 kpc beyond the break. They find that about 85% of 
these stars have formed in the inner disk and subsequently mi- 
grated outwards, beyond the truncation radius. The break in 
the surface density coincides with the minimum mean stellar 
age in the disk. Inwards of the break younger stars populate 
the larger radii, but the trend is reversed once we cross the 
point of truncation. In this scenario, the outermost regions of 
disks are old because the old stars have had the most time to 
migra te farthest ou t in the disk. While the positive age gradi- 
ent in|^^^Ferari()2008^ models could succee d in explaining 
the lGalleti e t al. (200i) and lBarker et al.1 ( 12007b results for the 
outer disk of M33, it is not clear whether the same mechanism 
would produce the old outer disk in the purely exponential 
disk of NGC 300. 

The question that arises here is whether what we observe 
as an exte nded stellar disk is in fact a halo component in 
NGC 300. Tikhonov et al] 12005') use HST/ACS to study the 
outer disk of NGC 300 and attribute the red giants found in 
their outermost field at 13' to a halo population. However, 
our surface brightness profile of NGC 300 in Figure [TO] con- 
firms that the outer disk indeed follows the same exponential 
profile as the inner disk and hence is not a distinct popula- 
tion. The possi b le rea sons for a discrepant result reported by 
iTikhonov et al.l (12005 ) are two-fold. Their outer field spans 
only 2' with large variations of star counts within the field, 
and they do not account for the contribution to the star counts 
by faint background galaxies. Correcting for the contamina- 
tion from the unresolved galaxy population would have the 
effect of steepening the profile in the outer field and making 
its slope more similar to the slope derived from their inner- 
most fields. 



Using the color of the stars on the red giant branch as 
a proxy for their metallicity, we derive the metallicity dis- 
tribution function of the stars in the outskirts of NGC 300 
and find that metallicity gradient changes slope at a radius 
of 15.5' 9.3 kpc). Does this mark the end of the disk 
in NGC 300 and a transition to an old halo? The argument 
against halo domination at large galactocentric radii is the rel- 
atively high metallicity in the outermost 10'. Typical stellar 
halos are more metal-poor than what we find for the outer 
disk of NGC 300. In the Lo cal Group, the halos of the Galaxy 
dBeers et al. 2005i ICarpUo et al. 2007; Ivezic et al. 2008) and 
M31 dChapman et al. 2006; Koch et al. 2007) have metaUic- 
ities of the order of [Fe/H] = -1.5, or lower. Even the less 
massive, almost bulge-less M33, believed to be a local ana- 
logue of NGC 300, ha s a metal-poor halo with [Fe/H] = -1 .3 
to -1.5 (iBrooks et al. 2004; McConnachie et al. 2006), and 
a more metal- rich disk with the photom etric metallicity of 
[Fe/H]w -0.9 (i McConnachie et al.ll2006l) . similar to the outer 
disk metallicities we derive. We also note here that stellar 
halos are usually found to exhibit a negative metallicity gra- 
dient (Harris 1996; Parmentieret al. 2000; Koch et al. 2007), 
which is in disagreement with our positive or flat gradient in 
the outer disk. However, due to the limited radial extent of 
our data, this is not as strong of an argument against the halo 
domination at large radii as are the high metallicities found in 
our outermost field. 

Negative abundance gradients are a common feature of 
disk galaxies. The dependence of yield, star formation and 
gas infall on the galactocentric rad ius in chemical evolution 
models reproduces the gradients (iGoetz & Koeppeiil 119921 : 
iMatteu cci & Francois 1989), which are shallow and with a 
large but uniform dispersion. A metallicity trend in which the 
decline in abundances with radius is followed by a metallicity 
plateau, and which we find in the outer disk of NGC 300, 
may be a general property of spiral disks. It has so far 
been observed in all of the galaxies for which the deep stel- 
lar photometr y is av ailable - M31, M33 and the Galaxy. 
IWorthev et all (l2005l) use archival HST images of a num- 
ber of fields spanning the galactocentric distances out to 50 
kpc to derive abundance gradient in the disk of M3 1 . While 
the inner disk exhibits a negative metallicity gradient, as ex- 
pected from the chemical evolution models, the gradient dis- 
appears at ~ 25 kpc with the metallicity re maining nea r ly cori - 
stant in the outer disk at \M/U] » -0 . 5. lYong et al.l (l2005l) . 
ICarnev et al] (|2005|) and lYong et al.l (l2006l) find a similar 
metallicity plateau at the radius of 10- 12 kpc in the outer disk 
of the Galaxy. They also observe the a-enhanced out er disk 
with [ a/Fe1 « 0.2. Similar results are those of T warog et alJ 
who find the break in the metallicity gradient of open 
clusters at R« 10 kpc, with the constant metallicity at larger 
galactocent ric distances . In the ir deep ACS study of the M33 
outer disk. Barker et al.l (12007) see an indication for the flat- 
tening in the metallicity gradient beyond the radius of 50'. 
The same effect has also been recently detect ed through the 
oxyge n abundances of H II regions in M83 (iBresolin et al.l 
l20Cil . 

We present here two scenarios to explain the metalhcity 
plateau (or upturn) in the outer disks of spirals. 

In the radial mixing scen ario, the mechanism described by 
ISellwood & BinnevI (|2002[) is responsible for transporting the 
stars within the disk while preserving t he nearly circular or- 
bits an d exponential light profile. In Bland-Hawthorn et alJ 
( 120051) we show that the outer disk of NGC 300 is a high-Q 
environment (with the Toomre Q parameter of ^ 5 ± 2) and 
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Fig. 14. — The illustration of a working model for the evolution of abun- 
dance gradients in disks. (1) Central spheroid is formed at high redshift 
and intense star formation estabhshes the steep metallicity gradient. (2) The 
metal-poor and uniform-metallicity disk is assembled at z ~ 2. (3) Star for- 
mation migrates to progressively larger radii and enriches the initially low- 
metallicity regions. (4) As a consequence, the abundance gradient continues 
to flatten. 

as such is expected not to be responsive to gravitational dis- 
turbances, such as star formation or spiral waves. If however 
there is a mechanism that could lower the estimated value of 
Q (e.g. presence of cold gas that is unaccounted for in our 
earlier calculation of Q which estimates the surface density 
from the H I observations), scattering of stars by spiral waves 
dSellwood & Preto..2002) could possibly explain the existence 
of stars in these diffuse regions. Potentially, stellar migra- 
tions could also explain the flattening or u pturn we observe 
in the metallicity gradient of NGC 300. Ro skar et al] (l2008h 
reproduce a similar change in trend of the stellar age with ra- 
dius. However, scattering over such large distances (5-8 kpc) 
would require spiral waves m uch stronger than those consid- 
ered by lSellwood & Binneyl (12002) . 

The accretion scenario is illustrated in Figure [141 In 
this picture, the central spheroid is assembled first in the 
early Universe (z > 4). The first generations of stars that 
form establish a steep abundance gradient in the spheroid 
(1), reaching the metallicities of around solar. This is sup- 
ported by solar metallicit i es see n in QSOs at high red- 
shifts (Hamann & F erlaiidl 1 19991 and references therein) 
and steep abundance gra dient i r i the old stellar bulge of 
the Galaxy JZo ccali et a ll 120081: iFeltzing & Gilmore 200Ct 
R-ogel et al."1999; Minni ti et all 1 199 51 butsee Ramire z et al] 
2000; Rich et al. 2007). The disk component is established at 
z ^ 2 t Wolfe et al.ll2005l) . The initial disk is metal-poor and 



with uniform abundance (2), reflecting the metallicities in the 
IGM at the formation redshift. As the galaxy is built in the 
inside-out fashion, the extent of star formation migrates to- 
wards the outer disk and the parts of the disk that previously 
had pristine abundances are now forming metals (3). This is 
a consequence of the gas and stellar density building up in 
the outer disk. As a result, regions of the disk at progres- 
sively larger radii become enriched and abundance gradient 
flattens with time (4). Observations of a varie ty of chemi- 
cal gradients in PNe (Maciel et al. 20031 l2006h . open clus- 
ters lMeretanr2002; Chen et al. 2003), cepheids, OB asso- 
ciations (Maciel et al. 2005; Daflon & Cunha 2004) and H 11 
regions (.Maciel et al.,,2007 ) have shown that older objects ex- 
hibit steeper abundance gradients resulting in chemical gra- 
dients that become flatter with time. We note however that 
while some of the chemical evolution models that incorpo- 
rate inside-out g rowth of galactic disks are able to reproduce 
this behavior (iMolla et al] 119971; [Bp issier & PrantzosI [19991; 
iPortinari & Chiosil Il999t ]Hou et al . 2000), others predict 
steep ening of gradient with time ( ^siill988i; Chiappin i et alJ 
1200 lb . This demonstrates the sensitivity of abundance gra- 
dients on adopted form of the radial dependence for the star 
formation and gas infall. 

A simple consequence of this picture is that there ex- 
ists a radius in a galactic disk at which a transition occurs 
between the star forming disk with a negative abundance 
gradient and a uniform-abundance initial disk, as observed 
here in the case of NGC 300, as well as in the disks of 
the G alaxy ( Yong et al. 2005; Carn gy et alJl2005t lYong et alJ 
120061) and M31 (.Worthey et alJi200 5'). The outer disk abun- 
dances reflect the local metallicity floor, which is supported 
by enhanced [ct/ Fe] ratios and low [Fe/H] both in exter- 
nal gas locaUy (IColUns et all I2003L l2007h and in the in- 
tergalactic medium at high redshift jProchaska et al.l 120031 ; 

mm4. f2me). 



iDessauges^ Zavadskv et alllool 12006') 

The accretion scenario suggests that stars today were born 
roughly in situ, whereas the radial mixing scenario implies 
that stars largely migrate away from where they were formed. 
Both of these have consequences for future studies of disk 
formation. 
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